Stripe and bubble phases in quantum Hall 

systems* 

Michael M. Fogler 
February 1, 2008 



Department of Physics, Massachusetts Institute of Technology, 77 Massachusetts 
Avenue, Cambridge, MA 02139, USA 



Abstract 

We present a brief survey of the charge density wave phases of a two- 
dimensional electron liquid in moderate to weak magnetic fields where 
several higher Landau levels are occupied. The review follows the chrono- 
logical development of this new and emerging field: from the ideas that led 
to the original theoretical prediction of the novel ground states, to their 
dramatic experimental discovery, to the currently pursued directions and 
open questions. 



1 Historical background 

Until recently, the quantum Hall effect research effort has been focused on the 
case of very high magnetic fields where electrons occupy only the lowest and per- 
haps, also the first excited Landau levels (LL). Investigation of the weak mag- 
netic field regime where higher LLs are populated, was not considered a pressing 
matter because no particularly interesting quantum features could be discerned 
in the magnetotransport data. The experimental situation has changed around 
1992, when extremely high purity two-dimensional (2D) electron systems be- 
came available. At low temperatures, T < 30mK, such samples would routinely 
demonstrate very deep resistance minima at integral filling fractions down to 
magnetic fields of the order of a tenth of a Tesla jl]. This indicated that the 
quantum Hall effect could persist up to very large LL indices, such as ~ 100, 
and called for the theoretical treatment of the high LL problem. Although 
the integral quantum Hall effect could be explained without invoking electron- 
electron interaction, two other experimental findings strongly suggested that the 
interaction is important at large N. One was the prominent enhancement of the 
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bare electron ^-factor g and the other was a pseudogap in the tunnehng den- 
sity of states Q . Surprisingly, no interaction-induced fractional quantum Hall 
effect has ever been observed at higher TV in contrast to the case of the lowest 
and the first excited LLs {N = and 1). An effort to understand this puzzling 
set of facts led A. A. Koulakov, B. I. Shklovskii, and the present author to the 
theory of charge density wave phases in partially filled N >2 LLs [Q, ^ . When 
this theory received a dramatic experimental support Q , a broad interest to 
the high LL physics has emerged. Below we give a brief review of this new and 
exciting field. Some of the ideas presented here are published for the first time. 
For previous short reviews on the subject see Refs. ||, M. 



2 Landau quantization in weak magnetic fields 

Consider a 2D electron system with the areal density rt in the presence of a 
transverse magnetic field B. For the case of Coulomb interaction, at zero tem- 
perature and without disorder, the properties of such a system are determined 
by exactly three dimensionless parameters: Ts, v, and Ez/HtUc- The first of 
these, Ts — (Trrta^)"^/^, measures the average particle distance in units of the 
effective Bohr radius as = fi^K/me^. The properties of the electron gas are 
very different at large and small . In these notes will focus exclusively on the 
case < 10. This is roughly the condition under which the electron gas in zero 
magnetic field behaves as a Fermi-liquid [Q. As we discuss below, the system 
is no longer a Fermi-liquid at any finite B; however, the basic structure of LLs 
separated by the gaps hwc, where Wc = eB/mc is the cyclotron frequency, sur- 
vives at arbitrary low B. In this situation, the second dimensionless parameter, 
u = 2TrPn, specifies how many LL subbands are occupied. Here I = {Hc/eBY^^ 
is the magnetic length. The lower LLs are fully occupied, while the topmost 
level is, in general, partially filled. Therefore, i' — 2N + vn, where the factor of 
two accounts for the spin degree of freedom and vn is the filling fraction of the 
topmost (A'^th) LL, < i/^ < 2 (see a cartoon in Fig. 

The remaining dimensionless parameter Ez/fnUc introduced above is the 
ratio of the Zeeman and the cyclotron energy. It affects primarily the dynamics 
of the spin degree of freedom, which is beyond the scope of these notes. Suffices 
to say that in the ground state the topmost A^th LL is thought to be fully 
spin-polarized for N > Jll| with a sizeable spin gap. In the important case of 
GaAs, the spin gap greatly exceeds Ez due to many-body effects. 

Because of the spin and the cyclotron gaps, the low-energy physics is domi- 
nated by the electrons residing in the single spin subband of a single (topmost) 
LL. All the other electrons play the role of an effective dielectric medium, which 
merely renormalizes the interaction among the "active" electrons of the A^th 
LL. This elegant physical picture was first put forward in an explicit form by 
Aleiner and Glazman [ p^ . 

The validity of such a picture in weak magnetic fields is certainly not ob- 
vious. Naively, it seems that as B and HlOc decrease, the LL structure should 
eventually be washed out by the electron-electron interaction. The following 
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reasoning shows that this does not occur (see also Ref. |Q for somewhat differ- 
ent arguments). Let us divide aU the interactions into three groups: (a) intra-LL 
interaction within iVth level, (b) interaction between the electrons of iVth level 
and its near neighbor LLs, and (c) interaction between A^th and remote LLs 
(with indices N' significantly different from TV). The last group of interactions 
is characterized by frequencies much larger than ujc- It is not sensitive to the 
presence of the magnetic field and leads only to Fermi-liquid renormalizations 
of the quasiparticle properties. Interactions within the groups (a) and (b) have 
roughly the same matrix elements but the latter are suppressed because of the 
cyclotron gap. Hence, it is the interactions among its own quasiparticles that 
are the most "dangerous" for the existence of a well-defined Nth LL. It is crucial 
that the intra-LL interaction energy scale does not exceed the typical value of 

E,^ ~ O.ley kR,, (1) 



where Rc — y/2N + 1 1 is the classical cyclotron radius IJ. The ratio 
Eex/hujc ''^ O.Irs is i?-independent; thus, there is a good reason to think that 
the validity domain of the proposed single-Landau-level approximation extends 
down to arbitrary small B's and, in fact, is roughly the same as that of the 
Fermi-liquid (r^ < 10). This is certainly borne out by all available magnetore- 
sistance data ||, ^, |^ . Henceforth we focus exclusively on the quasiparticles 
residing at the topmost LL. 

The inequality iJ^x < means that the cyclotron motion is the fastest 
motion in the problem, and so on the timescale at which the ground-state cor- 
relations are established, quasiparticles behave as clouds of charge smeared along 
their respective cyclotron orbits, see Fig. [^b. The only low-energy degrees of 
freedom are associated with the guiding centers of such orbits. In the ground 
state they must be correlated in such a way that the interaction energy is the 
lowest. This prompts a quasiclassical analogy between the partially filled LL 
and a gas of interacting "rings" with radius Rc and the areal density z^Ar/(27rP). 
Note that for i^jv > 1/-^ the rings overlap strongly in the real space. 

Strictly speaking, the guiding center can not be localized a single point, and 
so our analogy is not precise. However, the quantum uncertainty in its position 
is of the order of At large N, where I <C Rc, the proposed analogy becomes 
accurate and useful. For example, it immediately clarifies the physical meaning 
of E'ex as a characteristic interaction energy of two overlapping rings. 

Technically, the high LL problem is equivalent to the more studied N = 
case if the bare Coulomb interaction vo{q) is replaced by the renormalized 
interaction 



v{q) 



L 



N 



(2) 



where e(g) is the dielectric constant due to the screening by other LLs ^ 
and the bracketed expression compensates for the difference in the form-factor 

F^{q)=LJ^^)e-^'^'/' (3) 
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Figure 1: (a) Landau levels. Darkened ellipses symbolize electrons 
and arrows — their spins (b) A quasiparticle at the iVth LL viewed 
as a ring-shaped object immersed into a medium formed by the 
filled lower LLs. 



of the cyclotron orbit at iVth and at the lowest LLs, with Ln{z) being the 
Laguerre polynomial. In the next section we will discuss the consequences of 
having such an unusual interaction. 

3 Charge density wave instability 

The mean-field treatment of a partially filled LL amounts to the Hartree-Fock 
approximation, first examined in the present context by Fukuyama et al. It 
is worth pointing out the differences between their paper and our own 
workjj] reviewed below in this section. The pioneering work of Fukuyama et 
al. |14| appeared in 1979. A few years later, after the discoveries of integral and 
fractional quantum Hall effects, it became clear that for the exception of dilute 
limit, the Hartree-Fock approximation is manifestly incorrect for the lowest LL 
case . By 1995 when we started to work on the high LL problem, Ref. ||lj] has 
been effectively shelved away. In contrast to Ref. |lj] , who did not try to assess 
the validity of the Hartree-Fock approximation, our theory of a partially filled 
high LL Q was based on this kind of approximation because it is the correct 
tool for the job. This point is elaborated further in Sec. ^. Another important 
difference from Ref. [Q is a parametric dependence of the wavevector g^, of the 
CDW instability on the magnetic field: we find cx instead of cx B~'^/'^ 
in the theory of Fukuyama et al. Finally, the physical picture of ring-shaped 
quasiparticles that guided our intuition is quite novel and applies only for high 
LLs. 
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Figure 2: Direct, exchange, and the total Hartree-Fock potentials 
in g-space for N = 5 and rs — 0.5. [Reproduced with permission 
from Fig. 2 of Ref. g (a)]. 



Within the Hartree-Fock approximation, the free energy of the system is 
given by 

Fhf = ^^UHF(9)|(A(<z))P+fcBT^(nxlnnx + (l-nx)ln(l-nx)), (4) 

q X 

where A{q) = {'^'^^'^ / LxLy)Y.x ^'"'"^ °^x+qyi^ /2^x~qyi^ /2 are guiding center 
density operators, and Ly are system dimensions, iflx) are creation 
(annihilation) operators of Landau basis states and nx = a}-^ax are their 
occupation numbers subject to the constraint ('^^) ~ 

LyVN/{2TTp). Here 

and below we assume that < i^at < 1 because the states with 1 < i^at < 2 
are the particle-hole transforms of the states with 2 ~ vm and do not require 
a special consideration. The Hartree-Fock interaction potential is defined by 
WHF(g) = (1 - 5qfi)uH{q) - upiq), where 

UH{q) = i{q)e-'^"''^^ - 2nPuH{qf) (5) 

are its direct and exchange components (tildes denote Fourier transforms) Q. 
Their g-dependence is illustrated in Fig. |^. 

The inherent feature of uhf{q) is a global minimum at 

q, « 2A/R,, (6) 

where it is negative, 'iiHF(9*) < 0. Within the Hartree-Fock theory, it leads to 
a charge density wave (CDW) formation at low enough temperatures T [Q. 
Indeed, at high temperatures the entropic term dominates and the equilibrium 
state is a uniform uncorrelated liquid with (nx) — i^N and (A(g)) = dgfiU^. 
At lower T, it is more advantageous to forfeit some entropy but gain some 
interaction energy by creating a guiding center density modulation (A^) ^ 
with wavevector q = q^. 
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Since the quasiparticles are extended ring-shaped objects, the actual quasi- 
particle density modulation in any of our states is given by the product of 
the amplitude A(g) of the guiding center density wave and the cyclotron orbit 
form-factor: 

P{q) = Hq)FN{q). (7) 

The physical electric charge modulation in the system is further suppressed by 
the additional factor of e{q) due to the screening by the lower LLs. A peculiarity 
of the 1 case is that is very close to the first zero go of UH{q), which 

it inherits from the form- factor Fpf{q). |^ On the one hand, this means that 
the physical electric charge modulation is always rather small — a few percent 
in realistic experimental conditions. On the other hand, it explains why this 
instability develops in the first place. Indeed, usually the direct electrostatic 
interaction is repulsive and dominates over a weak attraction due to exchange. 
In our system uniq) vanishes at the "magic" wavevector go because no charge 
density is induced, p{qa) =0. As a result, the exchange part dominates and 
gives rise to a range of g's around qo where the net effective interaction uhf(?) 
is attractive, which leads to the instability. 

The nodes of FN{q) responsible for the vanishing of p{q) exist for a purely 
geometric reason that the quasiparticle orbitals are extended objects of a specific 
ring-like shape. The size of the orbitals is uniquely defined by the total density 
and the total filling factor. These two facts make the position of the global 
minimum q^ very insensitive to approximations contained in the Hartree-Fock 
approach as well as many microscopic details, e.g., the functional form of e{q), 
thickness of the 2D layer, which affects vo{q), etc. At asymptotically large N 
where the rings are very narrow, Fj^{q) is closely approximated by a Bessel 
function Jo{qRc); hence, Eq. (H). Surprisingly, Eq. is quite accurate even at 
A^ - 1. 

The mean-field transition temperature T™^ can be estimated as the point 
where the stability criterion l/etot(?) < 1 of the uniform liquid state is first 
violated. Here etot(g) is the total dielectric function, including both the lower 
and the topmost LLs. A simple derivation within the (time-dependent) 
Hartree-Fock approximation gives 



etot(g) = e(g) < l + UH{q) 



2TTVkBT 



vn{1 - vn) 



(8) 



The T-dependence of this function a,t q — q^ is sketched in Fig. ||. 
The stability criterion leads to the estimate 

" " y^iiFiSl*)^ (9) 

For Ts ~ 1, A 3> 1, and vjsi = 1/2, it can be approximated by 

ksT^^iiyN = 1/2) ~ omhwc. (10) 

^In contrast, in the N = case studied by Fukuyama et al., uh(9) does not 
have nodes. 



6 



1 














Figure 3: T-dependence of the dielectric function of the Uquid 
state. The unstable region is hatched. 



Compared to the expression originally given in Ref. |Q(b), a certain 1/iV term 
is omitted here, out of precaution that the Hartree-Fock approximation, which 
is valid in the large- TV limit, may not have enough accuracy to describe 1/N- 
corrections. This caveat should always be kept in mind when using numerical 
estimates of ksT^^ , such as those reported in Ref. For the typical experi- 
mental situation §, Q], Eq. @ gives T^^ - IK. 

4 Mean-field phase diagram 
4.1 Charge density wave transition 

A more systematic way to study the CDW transition is via a Landau expansion 
of the free energy F in powers of the order parameters (A(q)) where q are 
restricted to the locus of the soft modes, |q| = q^,: 

CO n 

^ = E««(^) E Ui^M)- (11) 

n—2 qi+q2 + ...+qn— 4— 1 

Note that in the quasiclassical large-iV limit, the order parameter (A(r)) is 
proportional to the local filling factor, (A(r)) — i/Ar(r)/27rP. 

The above linear stability criterion is equivalent to the condition 02 > 0. At 
= \ where the cubic term vanishes by symmetry, aj^ = 0, the Landau theory's 
estimate for coincides with Eq. (^) . It predicts a second-order transition |^ , 
|l8| , which occurs by a condensation of a single pair of harmonics, whose direction 
is chosen spontaneously, e.g., q — ±(7*x. The resultant low-temperature state 
is a unidirectional CDW or the stripe phase. Away from the half-filling, 03 7^ 0. 
Hence, at vm ^ i the transition is of the first order, takes place at a temperature 
somewhat higher than predicted by Eq. (^), and is from the uniform liquid into 
a CDW phase with the triangular lattice symmetry the bubble phase, 

see Fig. I (left). 

Near its onset, the CDW order brings about only a small modulation of the 
local filling factor, so that the topmost LL remains partially filled everywhere. 
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Figure 4: Left: Mean-field pliase diagram. Right: Guiding cen- 
ter density domain patterns at T = 0. Shaded and blank areas 
symbolize filled and empty regions, respectively. 



As T decreases, the amplitude of the guiding center density modulation increases 
and eventually forces expulsion of regions with partial LL occupation. The 
system becomes divided into (i) depletion regions where (A(r)) — and local 
filling fraction is equal to 2iV, and (ii) fully occupied areas where (A(r)) = 
(27r/^)~^ and local filling fraction is equal to 2N + 1 (however, see Ref. Q for 
a discussion of truly small r^). At these low temperatures the bona fide stripe 
and bubble domain shapes become evident, see Fig. ^ (right). 

4.2 Stripe to bubble transition 

Near Tc, the Landau theory predicts the stripe-bubble transition to be of the 
first order. This seems to be the case at T = as well, at least at large iV, where 
the this transition occurs at vn ~ 0.39 Q, see Fig. In systems with only 
short-range interactions a density-driven first order transition is accompanied 
by a global phase separation. An example is the usual gas-liquid transition. The 
densities Us and rif, of the two co-existing phases are determined by Maxwell's 
tangent construction. Fig. ^d. 

In the present case the long-range Coulomb interaction changes the situation 
drastically. The macroscopic phase separation into two phases of different charge 
density is forbidden by an enormous Coulomb energy penalty. Only a phase sep- 
aration on a finite lengthscale, i.e., domain formation may occur. Since stripes 
and bubbles are two the most common domain shapes in nature ||2^ , it opens an 
intriguing possibility that bubble- or stripe-shaped domains of the stripe phase 
inside of the bubble phase may appear, i.e., the "superbubbles" (Fig. ||) or the 
"superstripes" . Their size would be determined by the competition between the 
Coulomb energy and the domain wall tension 7. The superbubbles would have 
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Figure 5: (a) The energy density e in units of hwc/lirP as a 
function of z/jv for the bubbles (thick solid line), the stripes (dashed 
line), and the uniform liquid (thin solid line) for rs = 1, ^ 1, and 
T = 0. (b) Schematics of the conventional tangent construction. 
As explained in the main text, it is too crude to account for the 
specifics of the long-range Coulomb interaction. For example, for 
the graph on the left, such a construction would give a vanishing 
rib and the maximum possible rig, which is misleading. 



a diameter 

1/2 



and increase the net energy density by ~ 7p/a, P being the fraction of the 
minority phase. It turns out that this additional energy cost shrinks the range 
of the phase co-existence considerably compared to that in the conventional 
tangent construction. It may totally preclude the phase co-existence if 

— >K--.)^^. (13) 

where /s(6)(ft) is the free energy density of stripes (bubbles). 

Incidentally, these considerations are also relevant for the main transition 
from the uniform state into the CDW one at T — Tc- We can think of the pri- 
mary stripe and bubble phases as examples of a frustrated phase separation. At 

= 1/2 the conventional tangent construction would predict fc^Tc = £'cx/4, 
whereas the actual (mean-field) transition temperature is lower [Eq. ([Tc|)] be- 
cause of the extra energy density associated with the edges of the stripes. 



Returning to the case of superbubbles, we estimate d f ^(jj-^ / dn ~ fs{b)l 



2 



../n^ and 7 ~ Ei^^R^ng dX T — Q. The criterion ( p^ becomes (n^ — nb)jng < 
Kc^/nl ^ \/N . Thus, we can be certain that the superstructures do not ap- 
pear at high LLs. The cases of small N or high temperatures require further 
study m. 
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Figure 6: A cartoon of the conjectured superbubble phase. 



4.3 Transitions caused by particle discreteness 

With the periodicity of the bubble phase set by the preferred wavevector , the 
area of unit cell of the bubble lattice is equal to So = 2-\/37r^/gJ. The number 
of particles per bubble is therefore M = SqVn /'^'n'l^ ■ Using Eq. (J^ we obtain [Q 

M « SmN, N:$>1. (14) 

It is natural to ask whether this formula should be taken literally, even when 
it predicts a nonintegral M. Strictly speaking, CDWs with a fractional number 
of particles per unit cell are not ruled out. However, we choose to ignore such 
an exotic possibility because early Hartree-Fock studies for the lowest Landau 



level 19 1 concluded that only the phases with integer- valued M are stable. In 
our own numerical studies only integral M were examined. The results are 
reproduced in Fig. |^ (see Ref. Q for details). We found that the ground state 
value of M never deviates from the prediction of Eq. (|l^) by more than unity 
at iV < 10 and all vn- This confirms the robustness of the optimal period g*. It 
also suggests a practical rule for calculating the optimal M: one should evaluate 
the right-hand side of Eq. (|lj) and round it to the nearest integer. 

Under the assumptions we made, M exhibits a step-like behavior with unit 
jumps at a finite set of filling fractions. They correspond to the first-order 
transitions between distinct bubble phases. Everything said earlier about the 
possible phase co-existence near the first-order transitions applies here as well. 
For example, we do not expect any "bubbles of bubbles" at large N although 
at moderate TV there is such a possibility. As decreases, the lattice constant 
of the bubble phase changes smoothly in between- and discontinuously at the 
transitions, but always remains close to in/^/^q,, « 3.0i?c- Only in the Wigner 
crystal state (M = 1) the lattice constant is no longer tied to the cyclotron radius 
but is equal to {iTr/y/SvNY^^l, which is much larger than Rc at <^ 1/N. 

One warning is in order here. While very useful, the Hartree-Fock approxi- 
mation and the Landau theory of the phase transitions do not properly account 
for thermal and quantum fluctuations in 2D. A preliminary attempt to include 
these effects reveals additional phases and phase transitions of different order. 
The revised phase structure will be discussed in Sec. |^ below. 
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Figure 7: Cohesive energy E'^^^ in a set of bubble phases with 
different number of particles per bubble M. is defined as 

the interaction energy per particle ii'int relative to the uncorrelated 
liquid [where E^int — —{vn l'^)Ec^. Calculation parameters: A'' = 
5, = -\/2, T — Q. The crosses represent the Laughlin liquid 
energies. (Reproduced with permission from Fig. 2 of Ref. g). 

5 Validity of the Hartree-Fock theory 

It has been mentioned above that the Hartree-Fock ground state of a partially 
filled LL is always a CDW. In particular, for iV = it is a Wigner crystal |p^ . 
It is well established by now that the latter prediction is in error at most filling 
fractions vjsi. Instead, Laughlin liquids and other fractional quantum Hall states 
appear [p^ . The Wigner crystal is realized only when vjq is very small (dilute 
particles) or very close to 1 (dilute holes). Below we present heuristic arguments 
and numerical calculations, which together make a convincing case that the 
situation at large TV is different, so that the ground state has a CDW order at 
all and not just in the dilute limit. 

5.1 Quantum Lindenmann criterion 

A heuristic criterion of the stability of periodic lattices is the smallness of the 
fluctuations about the lattice sites compared to the lattice constant. Unless Ts 
is extremely small, the stripes and bubbles are contiguous completely filled (or 
completely depleted) regions of the topmost LL. In this situation the fluctuating 
objects are the edges of the stripes and bubbles. Using the suitably modified 
theory of the quantum Hall edge states (see Sec. ^), one can estimate their 
fluctuations 5r to be of the order of the magnetic length, 6r ^ I. Since the 
period A of the stripe and bubble lattices is at least a few Rc = \/2N +11, the 
Lindenmann criterion ^ A is well satisfied, and so the CDW states should 
be stable at any and sufficiently large N . Basically, when the amplitude of 
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the local filling factor modulation is appreciable and the stripes (bubbles) are 
so wide, they are very "heavy," quasiclassical objects and quantum fluctuations 
are unable to induce a quantum melting of their long-range crystalline order. 
This is in contrast to the small N case where the quantum fluctuations in the 
real space are of the order of the lattice constant except in the dilute limit. 



5.2 Diagrammatic arguments 

In a work jl^ published soon after our original papers |^ , Moessner and Chalker 
systematically analyzed the perturbation theory series for the partially filled 
high LL problem. They were able to achieve definitive results under two simpli- 
fying assumptions: (a) there is no translational symmetry breaking and (b) the 
range R of the quasiparticle interaction is smaller than I. [This interaction is 
given by the inverse Fourier transform of vo{q)/e{q)]. Under such conditions the 
Hartree-Fock diagrams were shown to dominate in the N ^ 1 limit. This is a 
very important result. Even though it enables us to make controlled statements 
only about the high-temperature uniform state, it certainly enhances the cred- 
ibility of the Hartree-Fock results at all T. Due to the screening by the lower 
LLs embodied in the dielectric function e{q) [Eq. (0)], the effective range R of 



interaction turns out to be of the order of |12 gb ^ h^n/me^ — 1^21 vrl; thus, 
condition (b) is satisfied at N ^ r^T^, which is not too restrictive in practice. 

There is still an interesting albeit academic question of what happens at 
truly small and intermediate filling factors where 1 <C -/V ^ rj^, so that 
R^ I. It will be discussed shortly below. 



5.3 CDW vs Laughlin liquids 

Obvious competitors of the CDW states at simple odd-denominator fractions 
vn = l/(2fc -I- 1) are the Laughlin liquids. The interaction energy per particle 
in a Laughlin liquid can be found summing a rapidly converging series |^l|] 

i^LL = -^i?c. + ^f:CKyK, (15) 

where Vk are Haldane's pseudopotentials 

Vk = ^ jd\uu{q)FK{V2q), (16) 

ii^cx (briefiy introduced in Sec. |^ as a characteristic energy scale) is defined by 

E,^ = Uex(0)/(27rZ2), (17) 

and ck are coefficients calculable by the Monte-Carlo method |2^, |j. Numeri- 
cally, about a dozen terms in the series (^) are needed to get an accurate value 
of E'll at vjs! = i and = i. In the large- limit we can also derive an 
analytical estimate, guided by the asymptotic relation 

Vk ^ Ucx(2\/if/0//^ K > 1. (18) 
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It indicates that -Ell is determined by the behavior of Uex{q) at g ~ Since 
the exchange and the direct interaction potentials are hnked by the Fourier 
transform [cf. Eq. (||)], has the effective range of R/P in the g-space [Q. 
Thus, two cases have to be distinguished. 

1. i? <C I. — Physically, this corresponds to ^ max{l, r~^} (recall that the 
interaction range R is of the order of as)- Provided R ^ I, only first few terms 
in the series (^5|) are important, leading to the estimate E'll ~ Uex{l/l)/P. On 
the other hand, the interaction energy per particle in the CDW ground state, 
EcDW ^ —Ucx{qo)/l'^ ^ — -Ecx, is significantly lower, i.e., the CDW wins. In the 
practical case of ~ 1, the CDW should be lower in energy than the Laughlin 
liquids at A'^ > Nc where Nc is a small number. The numerical calculations 
reviewed below indicate that this "critical" number is A^c = 2. 

2. R:^ I. — This regime appears in the parameter window 1 ^ <C r^'^. 
It is mostly of academic interest because it can be realized only in very high 
density 2D systems where Ts ^ 1. Such systems are unavailable at present. 

The theoretical analysis proceeds as follows. It turns out that Eq. ( |T8| ) is 
correct with a relative accuracy 1/lnA^ even for moderate K. With the same 
accuracy we can replace all Vk's in Eq. ( |l5| ) by 2itEcx- Using the sum rule 
J2k = (1 — '^Tv^)/'^ nil' arrive at the asymptotic formula 

^LL ^ -Scx/2, (19) 

which is of the same order as -Ecdw. To compare the energies of the two states, 
we have to exercise some care. Skipping the derivation, which relies heavily on 
another sum rule ||l^, for the Hartree-Fock states 

EI(A(q))p = ^.^, (20) 

we quote only the final result, 

£^CDW — — £^cx/2. (21) 

It signifies that the Laughlin liquid and the CDW have the same energy with a 
relative accuracy of 1/lnA^. We may understand this surprising near equality 
as follows. Consider a system of particles interacting via a long-range two-body 
potential v{r), which is nearly constant up to a distance R and then gradually 
decays to zero at larger distances. The particles are presumed to be spread 
over a uniform substrate "of opposite charge" with which they interact via a 
potential —nQv{r). This fixes the average particle density to be no. It is easy to 
see that for any configuration of particles, which is uniform on the lenghscale of 
R, two potentials nearly cancel each other: (i) the total potential created at the 
location of a given particle by all the other particles of the system and (ii) the 
potential due to the substrate. The net potential is equal to S = —v{0) because 
the particle does not interacts with itself. Due to the pairwise nature of the 
interaction, the interaction energy per particle (including the interaction with 
the substrate) is one half of S, i.e., —v{0)/2. Now we just need to recall from 
Sec. E that the energy of the self-interaction in our system is Ecx to recognize 
Eqs. (19|) and (|2l|) as particular cases of this general relation. 
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The above discussion has several implications. First, it clarifies the physics 
behind the arguments of Moessner and Chalker |l^ that the CDW states are 
likely to face a strong competition from certain uniform states when the interac- 
tion is sufficiently long-range, I. Second, it leaves the nature of the ground 
state at 1 <C iV <C rj^ an open question at the moment. The CDW states seem 
to be favored in numerical calculations done for 2 < < 10 and both ^ 1 
and fs <C 1 (see below). Truly high N have not been investigated yet. 

One intriguing possibility is the emergence of novel phases where the CDW 
ordering is not static but dynamic. One particular example is a quantum ne- 
matic phase, which can be visualized as a "soup" of fluctuating stripes. Such 
phases are actively discussed both in the context of the quantum Hall effect [^2[ 
p3| , m |5[ |6[ and the high-temperature superconductivity ^ . We 

would like to reiterate that an experimental search for these exotic phases would 
require very special samples, e.g., with much lower than presently available. 

5.4 Numerical results 

Trial wavefunctions. — The analytical estimates for E'll and i?cDW derived 
above become accurate only at very large N. At moderate N the comparison 
of trial states has to be done numerically. The procedure of calculating the 
energies of Laughlin liquids has been outlined above. It relies on the mapping 
of the A^th LL problem onto a problem at the lowest LL with the modified 
interaction, Eq. (||). To compute the energy of a Hartree-Fock CDW state we 
use the same trick: the trial state is chosen from the Hilbert space of the lowest 
LL, but the interaction potential is appropriately modified. 

The first step is to define the wavefunction of a single bubble: 

vI/o{M= n (-^-^.)xexp(-E^), (22) 

l<i<j<A/ \ i=l / 

where zj = xj + iyj are complex coordinates of M quasiparticles that compose 
this bubble. A well-known property of the Vandermonde determinant indicates 
that ^'o is in fact, a Hartree-Fock state. It is easy to see also that \l/o defines 
the most compact arrangement of M quasiparticles allowed at the lowest LL, 
i.e., a circular droplet of a completely filled LL centered at the point x = y = 0. 

The trial state we studied is composed of bubbles arranged in triangular 
lattice. It can be obtained by replicating the bubble (|2^) and translating its 
multiple copies to the appropriate lattice sites. This has to be followed by the 
antisymmetrization with respect to particle exchanges among different bubbles. 
The resulting wavefunction does not have a simple explicit form. Fortunately, 
to calculate the Hartree-Fock energy we do not need the wavefunction but only 
the particle density, see Eq. (^) and Ref. [|l9j. An excellent approximation for 
the latter is simply the sum of the densities of the individual bubbles. Strictly 
speaking, it is not a fully self-consistent solution of the nonlinear Hartree-Fock 
equations because of a small nonorthogonality among the wavefunctions of dif- 
ferent bubbles. However, even the nearest bubbles are effectively so far away 
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from each other that the deviations from the self-consistency are extremely 
small. If desired, the degree of self-consistency can be further improved using 
an iterative procedure of the relaxation type with the described density distri- 
bution as the initial guess. We have done this kind of calculations ^ and found 
that the iterations lower the energy of the state by less than one part in 10^, 
which does not affect the comparison with the Laughlin liquid energy (known 
much less accurately). 

On the basis of such calculations, we concluded that the CDW becomes the 
ground state at :^Ar = | for A/^ > 2 and a,t = ^ for > 3. For the latter 
fraction and N — 2 the energies of the two trial states are so close that no 
definite conclusion could be made. Nevertheless, in practice the samples always 
contain some amount of disorder, which would favor the CDW state over the 
liquid state. Therefore, we established A^c = 2 as the "critical" LL index where 
the transition from the Laughlin liquids to CDW phases occurs. Nc turned out 
to be the same both for ^ 1 and ^ 1, and whether or not we included the 
effect of the finite-thickness of the 2D layer |^ . 

Since the fractional quantum Hall effect (FQHE) is traditionally associated 
with the Laughlin states while the CDW does not exhibit the FQHE, our results 
imply that the FQHE is restricted to the lowest and the first excited LLs, N = 
and A'^ = 1. It is in principle impossible to observe the FQHE at N >2, and to 
date no one has. We are therefore led to propose the global phase diagram of 
the 2D electron systems shown in Fig. |[ 

Exact diagonalization of small systems. — Strong evidence in favor of the 
CDW order in a partially filled N > 2 LL has been given by Rezayi, Haldane, 
and Yang who studied systems up to 14 electrons by means of the direct 
numerical diagonalization of the Hamiltonian. Crucial for their success was 
employing periodical boundary conditions along the x and y-directions (torus 
geometry) . This setup avoids imposing the defects into the CDW lattice, aligns 
the CDW in a specific direction (which facilitates its detection), and enables 
to deduce the number of particles per unit cell simply from the multiplicity of 
the ground state manifold. Rezayi, Haldane, and Yang found no evidence of 
incompressible FQHE states at 2 < A^ < 6 and — i, i, |, etc. Instead, they 
reported the ground state degeneracies and quasi-Bragg peaks in the structure 
factor fully consistent with the formation of the stripe phase near the half-filling 
and a bubble phase at ^ and j. The periodicity deduced from the quasi-Bragg 
peak positions agreed with the Hartree-Fock prediction (^) within a few percent. 

Density matrix renormalization group. — Shibata and Yoshioka studied 
N = 2 case using another powerful numerical technique, the density matrix 
renormalization group. Although not exact, it is presumed to be highly accurate 
both for the ground state energy and the ground state wavefunction. They were 
able to study larger systems, up to 18 electrons. Shibata and Yoshioka presented 
pair correlation functions unambiguously showing the stripe and bubble phases 
and pinpointed the transition point between them to be at ~ 0.38. 
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Figure 8: The global phase diagram of a 2D electron system in 
the axes {N, vm) for the practical case Ts < 1 (schematically). The 
labels "bl", "b2", etc. denote the bubble phases with 1, 2, etc. 
particles per unit cell. As the magnetic field decreases, the system 
traces a zigzag path through this diagram, which starts at the lower 
left corner and proceeds along continuous segments [N, 0) {N, 2) 
connected by discontinuous jumps {N,2) — + [N + 1,0). One such 
jump is shown by the dashed line with the arrow. In principle, 
points away from the zigzag path can also be sampled if the func- 
tional form of the bare interaction can be modified sufficiently 
strongly compared to the Coulomb law. 
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6 Experimental evidence for stripes and bubbles 



6.1 Resistance anisotropy 

The existence of the stripe phase as a physical reality was evidenced by a con- 
spicuous magnetoresistance anisotropy observed near half-integral fractions of 
high LLs ^. This anisotropy develops at low temperatures, T < O.IK, and 
only in very clean samples. The anisotropy is the largest a,t v = 9/2 {N — 2, 
i^N = 1/2) and decreases with increasing LL index. At T = 25 mK it remains 
discernible up to ~ 11^ whereupon it is washed out, presumably, due to resid- 
ual disorder and/or finite temperature. The main anisotropy axes seem to be 
always oriented along the crystallographic axes of the GaAs crystal: [110] (the 
high-resistance direction) and [110] (the low-resistance one). Especially striking 
are the data obtained using the square samples in the van-der-Pauw measure- 
ment geometry: at = 9/2 the resistances along the two principal directions 
differ by three orders of magnitude. However, one has to keep in mind that the 
van-der-Pauw measurements exaggerate the bare anisotropy of the resistivity 
tensor |33| , due to current channeling along the easy (low- resistance) direction. 
Indeed, Hall-bar measurements, which provide a faithful representation of the 
resistivity tensor, show much smaller but still a very significant anisotropy up to 
7:1. Once the current distribution effects are taken into account, the following 
picture emerges. At T > O.IK, the transport is isotropic. As T decreases, the 
resistivity along the [110] direction {pyy) decreases but only slightly. In contrast, 
the resistivity in the [110] direction (pxx) rapidly increases, growing by almost 
an order of magnitude when T drops down to 25 mK. The anisotropy is the 
largest at vn — 1/2 but persists in a sizable interval 0.4 < vm < 0.6 of filling 
factors. Thus, as a function of vn, pxx exhibits a peak resembling the quantum 
Hall transition peaks in dirtier samples at higher temperatures. 

In contrast, the magnetotransport measurements near half-integral fillings 
of TV = and = 1 LLs reveal no significant anisotropics and no peaks in the 
longitudinal resistance. Instead, the resistance exhibits a minimum as a func- 
tion of vn, which deepens at T decreases Q. Such unambiguous distinctions 
indicate that the electron ground states at high LLs N > 2 are qualitatively 
different from those in lower LLs. The = 9/2 is the fraction that demarcates 
the transition to the realm of novel high LL physics. 

The emergence of the anisotropy is very natural once we assume that the 
stripe phase forms. It is based on two concepts: pinning of stripes by disorder 
and the edge-state transport. Each of these topics deserves a separate discussion, 
which will be given (in a brief form) in Sees. |l^ and ^ respectively. Here we only 
sketch the basic ideas. The pinning serves the purpose of preventing the global 
sliding of the stripes. This singles out the edge-transport as the only viable 
mechanism of current propagation. The edges of the stripes can be visualized 
as metallic rivers, along which the transport is "easy." The charge transfer 
among different edges, i.e., across the stripes, requires quantum tunneling and 
is "hard" because the stripes are effectively far away. Thus, if the stripes are 
preferentially oriented along the [110] direction, the sample would exhibit the 
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anisotropy of the kind observed in the experiment. 

The physical mechanism responsible for the alignment of the stripes along 
the definite crystallographic direction is debated at present (see Sec. 12). In 
the absence of external aligning fields, the stripe orientation would be chosen 
spontaneously. On the other hand, due to enormous collective response of the 
stripe-ordered phase, the stripes can be easily oriented by a tiny bare anisotropy 
of the medium or the substrate. 



6.2 New insulating states 

The existence of the bubble phases at high LLs is supported by another striking 
experimental discovery: reentrant integral quantum Hall effect (IQHE) at i/ w 
4.25 and « 4.75. The Hall resistance at such filling factors is quantized at the 
value of the nearest IQHE plateau, while pxx and pyy show a deep minimum 
with an activated temperature dependence. The transport is isotropic in these 
novel insulating states, pxx ~ Pyy The current-voltage {I-V) characteristics 
exhibit pronounced nonlinearity, switching, and hysteresis. Such phenomena 
are hallmarks of the glassy behavior common for pinned crystalline lattices and 
conventional CDWs. Hence, these observations are fully consistent with the 
theoretical picture of a bubble lattice pinned by disorder. At TV = 2 we expect 
only two bubble phases: with two {M = 2) and with one {M = 1) particle 
per bubble. Both phases are subject to pinning and should be insulating at 
T = 0. To understand the reentrancy phenomenon we have to take into account 
the finite-temperature effects. The M = 2 bubble phase is more rigid than the 
M = 1 (Wigner crystal) phase, and remains stable at temperatures where the 
Wigner crystal is already melted. As vn is varied towards the nearest integer, 
the insulating M — 2 bubbles are replaced by a conducting plasma formed in 
place of the Wigner crystal. Close enough to the integer i'n, the plasma is so 
dilute that weakly interacting quasiparticles become localized by disorder and 
the conventional IQHE results. 

Very recently, several more reentrant insulating states has been discovered 
also in the iV = 1 LL Their nature remains to be determined but it is 

tantalizing to suggest that these are also the bubbles phases. 



6.3 Other experimental findings 

It was shown in a set of remarkable experiments that the anisotropy near half- 
integral fillings can be strongly affected by an in-plane magnetic field . When 
applied along the easy resistance direction, the hard and the easy anisotropy 
axes interchange at iJy > 0.5 T |3^, When applied along the hard direction, 
the infiuence of is much less pronounced and is to somewhat suppress the 
anisotropy. This intriguing behavior is thought to originate from the orbital 
effects of in a finite-thickness 2D layer 1 17, 38| . They can be crudely described 



as squeezing of the cyclotron orbits in the direction perpendicular to -B|| . For 
such distorted orbits the stripe phase energy depends on the orientation of 
the stripes with respect to the in-plane magnetic field. Calculations based on 
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a suitably generalized Hartree-Fock theory of the previous sections show that 
the preferred orientation of the stripes can be both parallel and perpendicular 
to the in-plane magnetic field, depending on microscopic details of the real 
systems Q . For the specific parameters believed to accurately describe the 
samples examined in Refs. [^6|, the perpendicular orientation is preferred (in 
agreement with the experiment). However, further theoretical and experimental 
work is needed to fully understand these issues. 

The magnetotransport anisotropy at high LLs was also observed for p-type 
GaAs samples j39| . 

The higher current transport regime near u — was investigated. Gradual 
increase in the differential resistance along the hard direction was reported. 
Compared to the strong nonlinearities at the reentrant IQHE states, it is a 
relatively weak effect. 

Finally, the degree of anisotropy and the effect of the in-plane fields were 
found to be more pronounced in the lower spin subtend of the same Landau level. 
A possible explanation within the Hartree-Fock theory was recently suggested 
by Wexler and Dorsey . 



7 Many faces of the stripe phase 

In the wake of the experiments, a considerable amount of work has been devoted 
to the stripe phase in recent years ||2|, |4|, |44[ ||, ||, ||, |4|l. It led to 
the understanding that the "stripes" may appear in several distinct forms: an 
anisotropic crystal, a smectic, a nematic, and an isotropic liquid (Fig. ^). These 
phases succeed each other in the order listed as the magnitude of cither quantum 
or thermal fluctuations increases. Thus, at small N [N — 2, 3) or close to Tj, 
the phase diagrams of Fig. ^ and Fig. || need modifications to incorporate 
some (if not all) of those phases. The general structure of the revised phase 
diagram for the quantum (T = 0) case was discussed in the important paper of 
Fradkin and Kivelson [p4[ . Pinpointing the new phase boundaries in terms of the 
conventional parameters Ts and u will require further analytical and numerical 
work. Once again, we wish to emphasize that at large N these additional phases 
have very narrow regions of existence, if any. Let us now give the definitions of 
these intriguing phases and discuss their basic properties. 

Stripe crystal. — This state may in principle be understood at the Hartree- 
Fock level. It was shown Q| that the initially proposed Hartree-Fock solution 
with smooth edges and a strictly ID periodicity Q is not the global energy min- 
imum. A further gain in energy is attained once the stripes acquire a periodic 
modulation in the longitudinal direction, in antiphase on each pair of neigh- 
boring stripes. The resultant phase breaks the translational symmetry in both 
spatial directions and thus is equivalent to a 2D crystal. The unit cell of such a 
crystal is very anisotropic, with the aspect ratio of the order of : 1. However, 
there is only a single particle per unit cell; thus, this state is an anisotropic 
Wigner crystal. It is presumably the true ground state of the system at suffi- 
ciently large A^ where the Hartree-Fock is deemed to be exact. 
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Figure 9: Sketches of possible stripe phases. 

Smectic state. — The usual definition of the smectic is a "liquid with the ID 
periodicity." A smectic is less ordered than a crystal because the translational 
symmetry is broken only in one spatial direction. The rotational symmetry is of 
course broken as well. The smectic stripe phase can be thought of as a descen- 
dant of the stripe crystal where the longitudinal modulations on neighboring 
stripes persist locally, but have no long-range antiphase order because of dy- 
namic phase slips. In other words, the neighboring stripes are unlocked ]44| . In 
the thermodynamic limit this kind of state is equivalent to a stripe phase with 
no modulations (smooth edges), akin to the original Hartree-Fock solution Q. 
The necessary condition for the smectic order is the continuity of the stripes. If 
the stripes are allowed to rupture, the dislocations are created. They destroy 
the ID positional order and convert the smectic into the nematic. 

Nematic state. — By definition, the nematic is an anisotropic liquid. There is 
no long-range positional order. As for the orientational order, it is long-range at 
T = and quasi-long-range (power-law correlations) at finite T. The nematic 
is riddled with dynamic dislocations. Other types of topological defects, the 
disclinations, may also be present but remain bound in pairs, much like vortices 
in the 2D X-Y model. 

Isotropic liquid. — Once the disclinations in the nematic unbind, all the 
spatial symmetries are restored. The resultant state is an isotropic liquid with 
short-range stripe correlations. As the fluctuations due to temperature or quan- 
tum mechanics increase further, it gradually crosses over to the "uncorrelated 
liquid" where even the local stripe order is obliterated. 
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8 Effective theories of the stripe phase 



As often the case, the low-frequency long-wavelength physics of the system is 
governed by an effective theory involving a relatively small number of dynamical 
variables. The basic form of the effective theory is essentially fixed by the 
symmetry considerations. Let us outline how such theories are constructed for 
the various stripe states introduced in the previous section. 

8.1 Stripe crystal 

The low-energy dynamical variables of this state are the elastic deformation 
u = {ux{r,t),Uy{r,t)} of the crystalline lattice and the effective description is 
basically the elasticity theory. In addition, we have to account for the long-range 
Coulomb interaction between the density fluctuations n(r, t) = — noVu, where 
no = i^jv/(27rP) is the average particle density at the A^th LL. The symmetry 
arguments identify four nonvanishing elastic moduli: cn, C12, C22, and C44, so 
that the effective Hamiltonian takes the form 

+ in2(Vu)uH(Vu), (23) 

where mh should be understood as the integral operator. The dynamics of the 
system is governed by the Lorentz force and can be studied with the help of the 
effective Lagrangean 

£ — mnouJcUydtUx — H. (24) 

Solving the corresponding equations of motion, we find the following excitation 
spectrum of lattice vibrations (magnetophonons): 

c.(q) ^ ^ 

Here ujp{q) = [noUE_iq)q'^ /m]^^^ is the plasma frequency and 9 = B,rcta,n{qy / q^) 
is the angle between the propagation direction and the x-axis. For Coulomb 
interactions <jJp{q) ex. ^Jq leading to the well-known dispersion relation pO| 
io{ci) (X g^/^. In a strongly anisotropic stripe crystal, 012,044 <C cn causing 
the angular dependence w(q) oc sin 20 starting from relatively small q. All 
these results are valid for an idealized clean system. In reality the low-q magne- 
tophonon modes will be profoundly affected by disorder, which will be discussed 
in Sec. |l^. 

8.2 Smectic state 

As mentioned above, the smectic state is most closely related to the original 
Hartree-Fock solution Q]. It may be visualized as Hartree-Fock stripes slightly 
decorrelated by phonon-like thermal fluctuations. 



C44 + (cii + C22 - 2ci2 - C44) sin 26* 
4mno 



(25) 
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Harmonic approximation. — In the smectic only Ux retains its direct mean- 
ing of the elastic displacement. On the other hand, the density fluctuations n 
become an independent degree of freedom. For example, in the case of incom- 
pressible stripes, n comes from the stripe width fluctuations, which are separate 
from the "shape" fluctuations described by Ux- The number of dynamical vari- 
ables in the smectic and in the crystal is therefore the same. Moreover, the 
smectic can be thought of as a crystal that lost its shear rigidity because of 
phase slips between nearby crystalline rows. This intuitive picture enables us 
to deduce the effective theory for the smectic from Eqs. ( p3|) and ( p^ ) by a 
certain reduction. Of course, at the end we should verify that we did not miss 
any terms allowed by symmetry. The first step is to formally reintroduce Uy 
as a solution of the equation dyUy = —n/no — dxUx and use it to replace all 
instances of dyUy in Eq. (p3|). To ensure that the stripes are free to slide with 
respect to each other in the y-direction, the terms that depend on dxUy should 
be dropped: ci2 0, C44 0. Yet we need to be careful and recall that 
the complete elasticity theory always contains higher gradients such as (dyUx)^ 



[suppressed in Eq. (23)]. Once the coefficient in front of the first-order gradient 
term {dyUxY vanishes, higher gradients become dominant and must be included. 
The resultant effective Hamiltonian and the Lagrangean take the form 

H = ^{dxuf + ^{dluf + in(uH + X'^)n + CndxU, (26) 

C—pdtU — H, dyp = —mujc{n + n()dxu). (27) 

Here we switched to notations more natural for the smectic: Ux became u, Uy 
was traded for the canonical momentum p, the sum cn -|- C22 became Y , etc. 
The physical meaning of new phcnomenological coefficients is as follows. Y and 
K are the compression and the bending elastic moduli, x is the compressibility, 
and C — 2nl^Y dh\ K / dv M accounts for the dependence of the mean interstripe 
separation A on the average filling factor. 

Since the number of dynamical variables in the smectic is the same as in 
the crystal state, the collective mode count is also unchanged. We will keep 
referring to them as magnetophonons. Solving the equations of motion for n 



and u we obtain the dispersion relation of such magnetophonons |45 : 

1/2 

(28) 



c.(q) - "^('^^ '^^ 



LUc q 



Yql + Kqt, 



Unless propagate nearly parallel to the stripes, Loici) is proportional to sin 26 q'^/'^. 
Unlike in the stripe crystal, this relation is obeyed even at g ^ (again, in the 
absence of disorder or orienting fields). One immediate consequence of this dis- 
persion is that the largest velocity of propagation for the magnetophonons with 
a given q is achieved when 9 = 45°. 

Thermal fluctuations and anharmonisms . — From Eq. ( p^ ) we can readily 
calculate the mean-square fluctuations of the stripe positions at finite T, e.g., 

(KO.O) - ui0.yff) . - e.'.., . J^U ,29, 
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Figure 10: Portraits of the stripe phase on different lengthscales. 



This formula is valid if y is large so that magnetophonons with wavevectors 
Qy ^ can be treated classically [huj{q) ^ /csT]. As one can see, at any 
finite temperature magnetophonon fluctuations are growing without a bound; 
hence, the positional order of a 2D smectic is totally destroyed jS^ at sufficiently 
large distances along the y-direction, |y| ^ A\/YK /ksT = ^y. Similarly, along 
the x-direction, the positional order is lost at lengthscales larger than = 

{Y/Kfi^ey 

Another type of excitations, which decorrelate the stripe positions are the 
aforementioned dislocations. The dislocations in a 2D smectic have a finite en- 
ergy Ed ~ K. At UbT <^ Ed the density of thermally excited dislocations is 
of the order of exp(— iSu/fesr) and the average distance between dislocations 
is Aexp{2kBT/ Ed) ■ At low temperatures ^x,S,y <C C_d; therefore, the 

following interesting situation emerges (Fig. |l^). On the lengthscales smaller 
than (or ^j,, whichever appropriate) the system behaves like a usual smectic 
where Eqs. (p6|-p8|) apply. On the lengthscales exceeding £^d it behaves^ like a 
nematic [^2| . In between the system is a smectic but with very unusual proper- 
ties. It is topologically ordered (no dislocations) but possesses enormous fluctu- 
ations. In these circumstances the harmonic elastic theory becomes inadequate 
and anharmonic terms must be included. The most important anharmonisms 
are captured in the following elastic Hamiltonian, which should be substituted 
in place of the first two terms in Eq. (Eq) |5l| : 



rr _Y 
rlpi — — 



dxu - ^(Vw)^ 



K 



(30) 



It can be easily checked that the expression inside the square brackets, which is 
the compressional strain, is invariant under rotations of the reference frame by 
an arbitrary angle (f). For example, if in the initial reference frame u — 0, then 
in the new frame u{x,y) = (1 — cos(/))x + sin (f>y so that dxU — i(Vu)^ = 0. 

What is the role of anharmonisms? As shown by Golubovic and Wang |^ , 
they cause power-law dependence of the parameters of the effective theory on 

■^In a more precise treatment |53|, the lengthscales ^jj^ oc and 5i5u oc 



S,'^'^ are introduced such that S,DxS.Dy = 
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the wavevcctor q: 

for < (x\qy^y?^^, Qy < (y \ and 



K^Koi^yqy)-'/^, (31) 

^^^o(^.g.)"'/', (32) 



for qx ^ ^ and <C ?j7^(<Z2;C2:)^^'^- The lengthscale dependence of the pa- 
rameters of the effective theory is a common feature of fluctuation-dominated 
phenomena. Other famous examples include the criticality near phase tran- 
sitions and in systems at their lower critical dimension, such as the 2D X-Y 
model. It should be mentioned that the lower critical dimension for the smectic 
order is d = 3 js^, so that the 2D smectic is below its lower critical dimension. 
This is the reason why the scaling behavior (|3l|) and (^2|) does not persist indef- 
initely but eventually breaks down above the lengthscale £,d where the crossover 
to the thermodynamic limit of the nematic behavior commences. 

Equations (|3^) and (^2|) indicate that the compression modulus decreases 
while the bending modulus increases as the lengthscale grows. This can be un- 
derstood from the following qualitative reasoning. Thermal fluctuations create 
a lot of wiggles on the stripes. When an external compressional stress is applied, 
it can be relieved not just by compression but by flattening of the stripes. Since 
the latter involves unbending of the crumpled stripes and the bending costs 
less energy than the compression (g^ instead of q^), the apparent compressional 
modulus Y is smaller than its bare value Yq. Similarly, when one attempts to 
bend crumpled stripes, some compression is necessarily involved, and so the 
bending modulus K appears larger. 

The scaling shows up not only in the static properties such as Y and K but 
also in the dynamics. The role of anharmonisms in the dynamics of conventional 
3D smectics has been investigated by Mazenko et al. and also by Kats and 
Lebedev jS^. For the quantum Hall stripes the analysis had to be done anew 
because here the dynamics is totally different. It is dominated by the Lorentz 
force rather than a viscous relaxation in the conventional smectics. This task 
was accomplished in Ref. The calculation was based on the Martin-Siggia- 
Rose formalism combined with the e-expansion below d — 3 dimensions. One set 
of results concerns the spectrum of the magnetophonon modes, which becomes 

^(q) ~ smecos'/'eii^qf/^'^^^^.f^. (33) 

ijJc^x V mno 

Compared to the predictions of the harmonic theory, Eq. (p8|), the g'^/^-dispersion 
changes to q^^"^. Also, the maximum propagation velocity is achieved for the 
angle 9 « 53° instead oi 9 = 45°. These modifications, which take place at 
long wavelengths, are mainly due to the renormalization of Y in the static limit 
and can be obtained by combining Eqs. (^8|) and (|3^). Less obvious dynamical 
effects peculiar to the quantum Hall smectics include a novel dynamical scal- 
ing of Y and if as a function of frequency and a specific g-dependence of the 
magnetophonon damping p5[ |. 
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Figure 11: Splay (a) and bend (b) distortions in a nematic. 

The latter issue touches on an important point. Our effective theory defined 
by Eqs. (|2^) and ( ]27| ) is based on the assumption that u and n are the only 
low-energy degrees of freedom. It is probably well justified at T ^ but be- 
comes incorrect at higher temperatures. The point of view taken in Ref. p5[ is 
that in the latter case thermally excited quasiparticles ("normal fluid") should 
appear and that they should bring dissipation into the dynamics of the magne- 
tophonons. 

Another intriguing possibility is for quasiparticles or other additional low- 
energy degrees of freedom to exist even at T = 0. Such more complicated 
smectic states are not ruled out and are interesting subjects for future study. 



8.3 Nematic state 

Much like loosing the shear rigidity due to phase slips converts a crystal to 
a smectic, loosing the compressional rigidity due to mobile dislocations can 
convert a 2D smectic into a nematic. The collective degree of freedom associated 
with the nematic ordering is the angle (j){r, t) between the local normal to the 
stripes N and the x-axis orientation. Due to inversion symmetry, N and — N, 
i.e., (j) and + tt are equivalent, that is why N is often referred to as the 
director rather than a vector The effective Hamiltonian for N is dictated 
by symmetry to be 

H^^^(VNr + ^\VxN\^. (34) 

The coefficients Ki and K3 are termed the splay and the bend Frank con- 
stants 51 . They control the cost of the two possible elementary types of direc- 
tor nonuniformity shown in Fig. |ll|. Note that in the smectic phase = —dyU. 
This entails the relation i^s ~ K between the parameters of the nematic and 
its parent smectic. On the other hand, the value of Ki is expected to be de- 
termined largely by the properties of the dislocations . The elastic part has 
a particularly simple form if A'l = X3, in which case Hi^j {Ky,/2){\7 cj))'^ just 
like in the X-Y model. 

Another obvious degree of freedom in the nematic are the density fluctua- 
tions n(r, t). A peculiar fact is that in the static limit n is totally decoupled from 
N, and so it does not enter Eq. (^). However, since the nematic is less ordered 
than even a smectic, the question about extra low-energy degrees of freedom or 
additional quasiparticles become especially relevant. We believe that different 
types of quantum Hall nematics are possible in nature. In the simplest case 
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scenario N and n are the only low-energy degrees of freedom. This type of state 
has been studied by Balents Q and recently by the present author . It was 
essentially postulated that the effective Largangean takes the form 

C=^j-\dtN)'-H. (35) 

(As hinted above, the full expression contains also couplings between i9(N and 
mass currents but they become vanishingly small in the long- wavelength limit). 
The collective excitations are charge-neutral fluctuations of the director. They 
have a linear dispersion, 



Lu{q) =q\jKijcos^e + K3-/sin^0, (36) 

and resemble spinwaves in the X-Y quantum rotor model. 

Also discussed in Ref. was a dislocation-mediated mechanism of the 
smectic-nematic transition within the framework of duality transformations de- 
veloped earlier by Toner and Nelson Toner and Fisher and Lee |]57| . 
This theory predicts the existence of a second excitation branch with a small 
gap. This gapped mode can be considered a descendant of the magnctophonon 
mode of the parent smectic. 

Very recently, Radzihovsky and Dorsey |^ formulated a different theory of 
the quantum Hall nematics, whose predictions disagree with our Eqs. ( |35| ) and 
(^6[). Logically, there are two possibilities. Either, as mentioned above, there 
are several distinct kinds of nematic states possible in nature or some of the 
theoretical constructions advanced in Refs. ^ are incorrect. To resolve 
these issues it is imperative to bring the discussion from the level of effective 
theory to the level of quantitative calculations. One promising direction is 
to investigate concrete trial wavefunctions of quantum nematics, e.g., the one 
proposed by Musaelian and Joynt |Q : 

* - n(^^- - - ^''•)' - «'] X exp(-^ |z,f /4?2). (37) 

j<k 3 

Here a is a complex parameter that determines the degree of orientational order 
and the direction of the stripes. This particular wavefunction corresponds to 
V — ^. (As explained earlier, it can also be used investigate the higher Landau 
level states with = _|)- Recently, the work in this direction was continued 
by Ciftja and Wexler |26[ . 

Other contributions to the theory of quantum Hall nematics have been fo- 
cused on a finite temperature case. They include a work of Fradkin et al pl| 
who investigated the role of external anisotropy field in the 2D X-Y formulation 
and also a paper by Wexler and Dorsey Q , where a quantitative Hartree-Fock 
analysis of the Toner-Nelson disclination unbinding scenario has been done. 

It is worth mentioning that in the quantum case the chain of transitions 
crystal — > smectic — > nematic isotropic liquid may or may not be realized in 
full. A finite-size study by Rezayi et al. |M suggests that the transition from 
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Figure 12: (a) Geometrical meaning of variables in the edge-state 
effective theory, (b) Backscattering. 



the smectic to an isotropic phase, as the interaction parameters are varied away 
from their N = 2 Coulomb values, can also occur via a first-order transition, 
without the intermediate nematic phase. The resultant isotropic state is highly 
correlated and has little in common with the "uncorrelated" Hartree-Fock liquid. 
The natural candidates for the isotropic state include a Fermi-liquid-like state 
of composite fermions and a Pfaffian (Moore- Read) state ||5^ . Understanding 
all these competing quantum orders and transitions between them remains a 
major intellectual challenge. In contrast, the phase structure at larger N, away 
from the "transition point" = 2, is much simpler and is adequately captured 
by the Hartree-Fock theory. 

9 Edge state models 

Another theoretical approach to the physics of the stripe phases is based on 
the edge-state formalism. The advantages of the edge-state models are two- 
fold. First, they allow to bridge the gap between the microscopic theory and 
the formulations based on the elasticity theory or hydrodynamics (albeit under 
some crucial simplifying assumptions). Second, they enable one to calculate not 
only collective but also single-particle properties, such as the tunneling density 
of states. On the other hand, being more specialized than the hydrodynamics, 
the edge models have a somewhat restricted domain of validity. The models 
considered so far |4^, ^ apply when the stripes are incompressible, i.e., 
when they are contiguous regions of vj^ = 1 separated by sharp boundaries 
from the regions with i^f^ = 0. In our opinion the edge state models are better 
suited for describing the crystal and the smectic states. In the nematic state 
stripes exhibit violent quantum (or thermal) fluctuations and the picture of well 
defined sharp edges is questionable. 

The first step in constructing the edge theory is to identify the low-energy 
degrees of freedom with the deviations of the stripe boundaries from their 
equilibrium positions — {j ± vn /2)K, see Fig. |l2|a. Such deviations induce 



27 



extra ID guiding center densities along the edges, = ±u'^/ {2t:P). The next 
step is to find an effective Hamiltonian H that governs the interaction among 
To this end it is convenient to introduce the Fourier transform 

PM) = A^exp(-zg,xi) / dye-'"" V.(y), (38) 

and define the center-of-mass stripe displacement Ux and the net local density 
fluctuation per unit area, 

Ma;(q) = 27r/2 [exp(|g^zyjvA)p+(q) + cxp(-|9^j/ArA)p_ (q)] , (39) 
7i(q) = [e^]){^qxVN Iy)p+{aL} - cxp(-§g^iyjvA)p_ (q)] . (40) 

The general symmetry requirements compel the long-wavelength part -ffo of the 
Hamiltonian to take the form (|2^), reproduced here for convenience: 

Y K 1 

Ho = jidxUxf + -^{dluxf + -n(uH + X^^)n + Cnd^u^. (41) 

The edge-state formalism enables one go further and obtain quantitative esti- 
mates for the phenomenological parameters Y , K, x, and C from the micro- 
scopic Hamiltonian. It should be clarified that the exact calculation of these 
parameters remains out of reach for the moment. Available derivations p3| , 
^ ^ ^ are certainly not exact. It can be shown that all of them 
rely on one specific approximation scheme, the time-dependent Hartree-Fock 
approximation (TDHFA). In the spirit of the discussion in the previous sec- 
tions, we expect the TDHFA to be quantitatively accurate at ^ 1 but only 
qualitatively correct in the case of current experimental interest, N ^ 2. Let us 
therefore proceed with the exposition of the field-theoretic edge-state formalism. 

Besides Hq, the full effective Hamiltonian contains another term i?2(jF : whose 
importance was first emphasized by Fradkin and Kivelson ||2^. To clarify its 
origin let us recall that in the Landau gauge A = (0, Bx, 0) each single-particle 
basis state \X) has a definite momentum in the y-direction, proportional to the 
mean value X of the x-coordinate in such a state, ky — X/P. In the simplest 
mean-field realization of the stripes, a given state is filled if ky e {x-'_^ / P , x-'_ / 
and is empty otherwise. Thus, ky — a;j_ play the role of Fermi points in an 
equivalent quasi-lD electron system. Note that "+" ("— ") stands for the sign 
of the Fermi velocity near a given Fermi points, in the positive (negative) y- 
direction. The density fiuctuations propagate along the edges only in the 
direction specified by this sign. Such a property is called chirality. Thus, the 
stripe phase is characterized by a Manhattan grid of alternating chiral edge 
states. From this perspective, Hq describes the particle-particle interactions, 
which cause only small changes in X and thus small momentum transfer in the 
vicinity of the Fermi points. Another low-energy process is to scatter particles 
between different Fermi points, e.g., x^^/P x''_/P and x-^^^/P x'^^/Z^, 
see Fig. Such scattering acts must involve a pair of particles to conserve 

the total momentum: one particle gains and the other looses the same amount 
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of momentum, 2qp = (x^ — x-'_)/P — vjssh.jl'^. H^q^ accounts precisely for the 
processes of this type. Since they cause the reversal of the propagation direction 
for the particles involved, it is natural to call it backscattering. Remarkably, it 
is possible to express ll2qF in terms of density fluctuations p^. Without 
going into details, we quote the final result 

- ^A"cos[(0i-0r') + ('^+-'/'-^')] (42) 
j 

Here are auxiliary dynamical variables related to as follows: 

p4 = 9,04/(2^). (43) 

The phenomenological parameters A° are proportional to ■UH(2g^) but in general 
depends on how the theory is formulated (the ultraviolet cutoff). In principle, 
interactions can scatter the particles not only between nearest edges but also 
next-nearest ones, etc. The corresponding amplitudes are proportional to mh(<z) 
where q is the appropriate momentum transfer. Since uh decreases exponentially 
at such large g, these other processes are expected to have negligible effect. Thus, 
the total Hamiltonian is H = Hq + H2qp ■ 

The final step in constructing the effective edge theory is determining the 



kinetic term. Either from Wen's bosonization theory |42 of by comparing 
Eqs. (27) and (Esh, one can come to the conclusion that (jA^ and are canoni- 



cally conjugate variables, so that the appropriate Lagrangean is 

^ = ^ E ^dtct^idy4>l -Ho- H2q, . (44) 

Let us now explain how this theory can lead either to crystalline or to smectic 
behavior. The idea is to treat H2qp as a small perturbation 
If this perturbation is irrelevant, in the long- wavelength low- frequency limit H 
reduces to Hq and C to the smectic form (p7|). On the other hand, if H2qp is rel- 
evant, then the and "— " edges of each stripe become strongly mixed so that 
a static 2qp density modulation in the y-direction appears. Its spatial period 
is a = 2t: / l2qp) — 27rZ^/z/ArA; hence, the number of particles on a given stripe 
within one modulation period is vn^ x a/{2TrP) ~ 1. This can be visualized as 
a ID chain of particles. The neighboring chains are locked in antiphase to lower 
their interaction energy. This state is an anisotropic Wigner crystal with the 
aspect ratio of the unit cell A : a = (i^Ar/27r)(A//)^ '--^ : 1. Two types of stripe 
crystals are possible. If the first term in H2qp is relevant, it is an "electron" 
crystal, if the second term is relevant, it is a "hole" crystal. 

The elasticity theory ( p3| ) of, e.g^ the electron crystal can be recovered as 
follows. Since the first sum in Eq. ( [l2|) is relevant, the large fluctuations of the 
argument of the corresponding cosines are inhibited. It is legitimate to expand 
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these cosines to arrive at 



- - E y("'. - - ^dy^^f, (45) 

j 

u,(q) EE -il.[</,+ (q)e'A9==/2 „ <^_(q)e-»A«^/2]. (46) 

In the long-wavelength limit it gives a term proportional to (dxUy+dyUx)'^ , which 
makes the shear modulus C44 finite [cf. Eq. (p3|)]. A more careful treatment 
presumably recovers the remaining elastic constant C12. Our choice of Uy in 
Eq. ( ^ ) is not arbitrary because we want to preserve the physical interpretation 
of Uy as the elastic displacement. This requires H2qp to be invariant under the 
shift by a lattice constant, Mj, — > Mj, + a. Since H2qp is invariant under the shift 
(/)^ — !■ (/)-!|_ + 27r, it fixes the coefficient of proportionality in Eq. (^6|) to the value 
given. 

The relevance or irrelevance of i?2gF depends on the stiffness of the stripes 
measured, roughly, by the dimensionless parameter yZ'*/'UH(7r/A). If it is small 
enough, then high-g and to fluctuations of and (either quantum or ther- 
mal) are sufficiently strong to cause effective averaging out of the cosines in 
H2qF ■ As a result, the coefficient A renormalizes to zero in the low q and lo limit. 
In contrast, for rigid stripes, the averaging of the cosine terms is not important. 
The perturbative renormalization group (RG) analysis ^ allows 

to make this argument more precise. However, to decide between the crystal 
or the smectic behavior, the RG requires accurate estimates for Y and other 
parameters of the edge theory as an input. Attempts to extract such parameters 
from the TDHFA for the case of current experimental interest, N ^ 2, led to 
contradictory statements in the literature In this regard, we wish to 

reiterate that at iV ^ 2 the stripe phase is very fragile and faces a strong com- 
petition from other quantum Hall states. Therefore, the Hartree-Fock and its 
time-dependent extensions are not quantitatively reliable. At the present stage 
the controlled theoretical analysis can be envisioned only for large N, building 
on the work of Moessner and Chalker Note that the ratio Y / U}i{Tr / A) 

becomes an A'^- independent number of the order of unity at A^ 3> 1 and ~ 1. 

Concluding this section, let us discuss the tunneling into the stripe phase 
from the normal metal. Two setups can be imagined: tunneling from a bulk 
metal or from an STM tip. In the first case the electron-electron interactions are 
screened by the metal, in the second they remain long-range. The differential 
tunneling conductance Gt — dl/dV at bias voltage V is proportional to the 
tunneling density of states g{eV), which is defined by g{E) = {1 / TT)QmG ii{E / H) , 
G[i{Lo) being the single-electron Green's function, 

GR{t)^-ie{t){^{0)^>\t)). (47) 



^ These four works are in mutual agreement on how the RG procedure needs 
to be formulated. A different type of RG analysis suggested in Ref. is 
believed to be in error. 
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Figure 13: The sketch of the true tunnehng density of states (sohd 
hne) and its Hartree-Fock approximation (dashed hue). 



Here and ^' are the electron creation and annihilation operators. The crude 
overall behavior of g{E) can be obtained at the Hartree-Fock level. In the sim- 
plest case of the smectic (unmodulated stripes), g{E) is determined by the slope 
of the self-consistent potential E{X) that defines the energies of the quasipar- 
ticle states \X), 

gHF{E) = \4TTfA{dE/dX)\-\ (48) 

The tunneling density of states calculated according to this formula is nonvan- 
ishing in a finite interval of width E^x around E — 0. At both ends of such 
an interval one finds two divergencies (van Hove peaks). The interior of the 
interval can be described as a shallow pseudogap, see Fig. |l3|. The negative- 
peak corresponds to tunneling into the states with X's in the middle of the 
filled stripes and the positive-i^ one — into the middle of the empty stripes. 
Conversely, at low energies (low V) the tunneling can only occur into points in 
a vicinity of stripe edges The Hartree-Fock results for g{E) are certainly 
not exact. However, the van Hove peaks and the pseudogap are expected to be 
true features and the Hartree-Fock estimate for the energy separation between 
the peaks should be quite reliable. The largest deviations from the Hartree- 
Fock predictions should occur at low E, which we now address. We will use the 
evolution in the imaginary time t — it picture common in quantum tunneling 
problems. 

Since the edges are effectively far away from one another, the incoming 
electron is initially accommodated into the closest edge, say, j+. It produces 
a compact charge disturbance with a high Coulomb energy. The charge has 
to gradually spread over the entire area. This takes place by the emission 
of the magnetophonons. (In the context of the edge-state theories they are 
traditionally called edge magnetoplasmons) . What is missing in the Hartree- 
Fock picture is precisely this relaxation from a point-like disturbance into an 
uniform state with an added quasiparticle of the given small energy E. Its net 
result is a suppression of the true g{E) compared to giip{E). To calculate the 
suppression factor, one can use the bosonization prescription '5 cx exp(i0'^) |^^ , 
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which entails 



GBXiiOn) oc / drexp itUnT - -(0^+(O)</)^+(t) - 0^+(O)0^+(O)) . (49) 



This expression can be evaluated with the help of the effective action ( p4| ) and 
then analytically continued to real frequencies iw„ — > uj+iO pv^. For the smectic, 
such a calculation has been performed by Lopatnikova et al. |49|| . They obtained 
I(V) oc exp(-VW^) and I(V) cx exp[- ln2(Vb/V^)] for the Coulomb and 
short-range interaction case, respectively. The last formula was also obtained 
by Barci et al. [Q. As one can see, the tunneling conductance vanishes at zero 
voltage, which is a typical result for tunneling into a correlated electron system. 
For the crystal, the calculation is more difficult because the action is essentially 
nonlinear [the harmonic approximation ( ^5|) is not adequate for the tunneling 
problem]. It has not been reported in the literature. On the physical grounds, 
we expect the tunneling conductance to be exactly zero unless eV exceeds the 
energy needed to nucleate a phase slip Uy Uy + a (an analog of an interstitial 
in conventional crystals). 



10 Pinning of stripes by disorder 

Disorder in the form of randomly positioned impurities is expected to strongly 
affect the degree of the positional, orientational, and topological orders in the 
stripe phase at low T. According to modern understanding, none of these sur- 
vives in the thermodynamic limit in two dimensions. However, in the most in- 
teresting case of weak disorder, the details of how each type of order is destroyed 
are subtly different. Let us start the discussion with the case of a smectic. Here 
the positional order is limited primarily by random elastic distortions. They 
can be characterized by the roughness function 

/(r)-(Kr')~w(r' + r)]2), (50) 

so that the positional correlation length can be defined by the condition 
/(^p) — h? . On the other hand, the correlation lengths for the topological 
order and orientational order are set by the mean distance between the 
disorder-generated dislocations and disclinations, respectively. Note that an 
appreciable magnetotransport anisotropy can be detected only if the sample 
size is smaller than ^o. This explains why samples of macroscopic dimensions 
L ^ 0.25 mm show no anisotropy unless they are extremely pure 0]. In 
addition to weakness of disorder, the gigantic in these experiments may have 
originated from a small bare anisotropy, which favored the alignment of the 
stripes along the [110] direction. Indeed, an external anisotropy qualitatively 
modifies the smectic state properties by generating a finite tilt modulus Y±. 
Correspondingly, the elastic part H^i of the Hamiltonian becomes 



Hp] — — 
2 



d^u + ^(V^)^]' + ^{dyuf + ^{dluf. (51) 
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The full Hamiltonian H = H^i + i?pm also includes the interaction with dis- 
order. It can be approximated by iJpm — J exp[iq^u{r)]V (r) where V{r) is a 
random short-range potential. The behavior of the static correlation functions, 
such as f{r), as functions of r and disorder strength has been evaluated within 
this model by Scheidl and von Oppen |^ building on the theory developed 
earlier for other pinned systems ^ . Scheidl and von Oppen suggested that 
in the current experiments the characteristic lengthscales satisfy the chain of 
inequalities < < i < Co, so that the stripe phase resembles an instanta- 
neous snapshot of a nematic with no long-range positional order and a number 
of frozen-in dislocations. 

Besides decorrellating the stripe positions, another interesting disorder ef- 
fect in the smectic phase can be inducing quenched random 2qi;'-modulations 
along the stripes, as suggested by Yi et al. This is similar to the effect of 
impurities on a ID "Wigner crystal" (where without impurities the density 
is uniform because of strong quantum fluctuation). To properly describe this 
effect and its consequences, a comprehensive quantum theory of pinned smectics 
is required. 

The pinning effects in the stripe crystal phase are qualitatively similar be- 
cause Hf,i is already not too different from the elastic Hamiltonian of an anisotropic 
crystal. The main difference is the doubling of the number of components of 
the elastic displacement field u — {ux,Uy\. 

The combined effects of finite temperature and disorder have not been in- 
vestigated in any detail, but based on the extensive work done in other con- 
texts ||6^, |6^ we may expect a thermal depinning at a certain "glass transition 
temperature." 

Exceedingly interesting topic is the dynamics of pinned electron (liquid) crys- 
tals in high magnetic field. First of all, there are number of open questions in the 
general theory of pinned manifolds (such as CDWs, vortex lattices, magnetic 
domains, etc.). The current interests in this subject area include dislocation 
dynamics and non-equilibrium states (e.g., high electric field sliding phases). In 
addition to that, there are interesting dynamical effects unique to quantum Hall 
systems. For instance, in the recent years it was experimentally discovered |l65| 
that even the linear response of a pinned Wigner crystal in a magnetic field is 
dramatically different from conventional expectations []66| . Typically, one antic- 
ipates the response to be dominated by low-frequency magnetophonon modes, 
which in the presence of pinning become localized wavepackets. Since the dis- 
order is different in different parts of the sample, a considerable inhomogeneous 
broadening is expected. Instead, a sharp resonance-like mode in the microwave 
absorption spectrum appears (typically around 1 GHz) . Assuming that the qua- 
siclassical description is adequate, Huse and the present author |^ explained 
the appearance of a narrow mode as a combined effect of the special Lorentz- 
force dynamics and long-range Coulomb interaction. However, the theory was 
done for the T — Q while the experimental line width seem to be dominated by 
thermal broadening. In fact, the thermal broadening is much more pronounced 
than one would naively expect and may reflect some nontrivial physics over- 
looked in a formulation based on harmonic elasticity theory p37|]. The similar 
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pinning resonance should exist in a stripe crystal as well and may even persists 
in a smectic phase. Since at high LLs all the energy scales are smaller, the 
pinning mode should have a lower frequency, perhaps, of the order of 100 MHz 
at = 9/2 for the kind of samples used in Refs. |^. 

11 Magnetotransport in the stripe phase 

The magnetotransport in a partially filled Landau level is a problem of two 
extremes. In a perfectly clean system, it is completely trivial and amounts to 
the uniform sliding motion of all the particles in the direction perpendicular 
to the applied electric field, leading to ctx^ — 0, cTxy — {e^ /h)v. On the other 
hand, in the presence of disorder, the magnetotransport immediately becomes 
very complicated. Especially for the stripe phase, the magnetotransport theory 
is very incomplete at the moment. For instance, the transport in the nematic 
phase is a total enigma. Certain progress has been achieved regarding the crystal 
and the smectic phases. As discussed above, pinning by random impurities 
eliminates the possibility of a global sliding motion of the electron liquid under 
the action of an external electric field. As a result, the stripe crystal is an 
insulator. In the smectic, gapless edge states may exist in which case the edge 
transport is possible. The quasiparticles travel along the edges in the direction 
prescribed by the edge chiralities and also be scattered between the edges by 
impurities. The scattering probability depends on the distance between the 
edges. For a general filling fraction the widths of the filled (electron) and 
empty (hole) stripes are unequal, and so the corresponding mean-free paths are 
also different, l^. ^ Ih- On a large scale the quasiparticle motion resembles a 
random walk with unequal elementary steps: along the stripes, it is a properly 
weighted average of 1^ and Ih, across the stripes it is the interstripe separation 
A. If we assume that (a) the stripes are continuous and are all oriented along 
the y-direction, (b) the scattering occurs only between nearest edges, and (c) 
quantum localization effects can be neglected, then a simple calculation leads 
to the following formulas for the components of the conductivity tensor |4|, ^ : 

A^i A / Z \ 

At the symmetric point le = lh^ the transport anisotropy ratios are given by 
^ — '^yyl'^xx = Pxxl Pyy = l\l ^ which is large for dilute impurities. However, 
there are reasons to doubt that Eq. ( |5^ ) directly applies to the experimental 
practice. Indeed, from the empirical estimate R ^ 7 one would conclude that 
le ~ 2.6A ~ 0.2 ^m, which seems too short for the extreme purity samples used. 
The most probable reason for reduction of the anisotropy ratio compared to 
the idealized formula (^2|), is the absence of high degree of the orientational 
and topological orders. As explained above, the pinning by disorder introduces 
pronounced elastic deformations as well as dislocations and disclinations into 
the perfect stripe pattern. Thus, in a more realistic model the above condition 
(a) should be abandoned. This leads to a view of the low-temperature stripe 
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phase is a collection of elongated puddles of j^jv = 1 state preferentially aligned 
in a certain direction (a kind of frozen nematic). The conductivity tensor can 
no longer be calculated in any simple way; however, it satisfies the generalized 
Dykhne-Ruzin semicircle law [ |69| 

Wyy + (fT.y - aof = {e^/2hf, ao = {e" /2h){2N + 1), (53) 

pointed out by von Oppen, Halperin, and Stern J68[. In addition, exactly at 



half- filling Gxy — cq and the product rule |69| ^ ^ applies: 

(JxxfTyy = {e^/2Kf. (54) 

Both relations agree quite well with the low-temperature experimental data at 
= 9/2 @|, which is encouraging. However, one has to be careful with drawing 



conclusions from such a comparison. The status of Eqs. ( pSf ) and (54) in the 
quantum domain is poorly understood and is even somewhat compromised by 
the lack of agreement with the data from dirtier samples. The situation is 
complicated by a unconventional behavior of the resistivity peak at = 9/2 
(low-T width saturation) , precisely where the product rule is supposed to apply. 

At higher temperatures where the CDW amplitude is small, the stripe edges 
are not well defined objects. It is better to describe the system as a "compress- 
ible liquid." A new type of transport theory is needed in this regime. It is 
reasonable to assume, for example, that the stripes can be depinned already by 
a vanishingly small electric field, so that a global sliding of the electron liquid 
becomes possible. If disorder is sufficiently smooth, it would be similar to a flow 
of a viscous 2D liquid across a disordered substrate jl6j . 

Sufficiently close to Tc the dislocations must become highly mobile so that 
the stripe phase should behave as a nematic on all relevant lenghtscales. A phe- 
nomenological approach to the magnetotransport in this regime was advocated 
by Fradkin et al. They argue on the basis of symmetry that the deviation 
of the resistivity ratio R = Pxxl Pyy from unity is proportional to the degree 
of the orientational order in the nematic. Under a simplifying assumption that 
the Frank constants are equal, they map the problem onto a 2D X-Y model in 
the presence of a small orienting field. Using the Monte-Carlo results for the 
latter model, they fitted the experimental data for R with the help of two free 
parameters. 



12 Other topics 

The space constraints force us to switch to a telegraphic-style overview of the 
remaining topics explored in connection with the CDW phases in high LLs. 

One of those is the bare anisotropy that aligns the stripes in the [110] di- 
rection. Its origin remains unclear. A careful study by Cooper et al. ff\ has 
ruled out a few earlier suggestions that included substrate morphology JTot , 
vicinal steps from a slight miscut of the wafer, and a small anisotropy of the 



effective mass 1 72 , 73] . One of the remaining viable alternatives is piezoelectric 
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effects IQ . The orienting effects of the in-plane magnetic field was further stud- 
ied in wide quantum wells [7^ . The combination of the empirical findings and 
the Hartree-Fock calculations suggest that the anisotropy is indeed tiny, ^ 1 mK 
per electron. The corresponding analysis for the p-type GaAs has also been 
performed | |75| . It is more involved because of the band-structure anisotropy. 
The preferred stripe orientation was predicted to sensitively depend on various 
details such as the width of the quantum well. 

The time-dependent Hartree-Fock calculations of the collective modes of the 
stripe crystal has been done by Fertig et al. |47). At low g's the results are in 
agreement with the general structure predicted by the effective theories, Sees. 
and p. At higher q a new feature is revealed: a roton-like minimum at g « g*, 
which may be thought of as a precursor to the stripe-bubble transition. 

Stripe phases in bilayer systems have been studied by Brey and Fertig |Q 
and also by Demler et al. . The additional layer degree of freedom leads to a 
rich phase diagram. Some of that physics may be operational in wide quantum 
wells @. 

The emergence of the CDW indicates some sort of instability of the con- 
ventional uniform FQHE states at high LLs. Such an instability was indeed 
found in numerical studies of Scarola et al. who relied on the composite- 
fermion theory. It was also discussed within the phenomenological mean-field 
composite-boson formulation . 



13 Conclusions 

The physics of high LLs has fledged into a fast growing and vibrant field. Despite 
a significant progress, a multitude of open problems remains. Some of them are 
brand new (the existence of qualitatively novel states of matter such as quan- 
tum nematic), others are venerable quantum Hall problems with a new twist 
(the integral quantum Hall transition in the stripe phase) . The real progress in 
understanding these intriguing yet difficult problems is likely to be tied with fu- 
ture experimental advances. It should be mentioned that the magnetotransport 
studies have certainly not exhausted their potential. Ongoing fantastic achieve- 
ments in sample fabrication continue to deliver the ever more perfect systems, 
which should expedite the discovery of new phases. A fresh example is the afore- 
mentioned group of new insulating states in the = 1 LL of 3.1 x lO'^cm^/Vs 
mobility sample ||3^ (Sec. Besides relying on increasing sample quality, it 
may be interesting to experiment with different sizes of already available sam- 
ples in search for the sequence of "spaghetti" phases sketched in Fig. ^ On the 
other hand, finite-frequency tools, such as microwaves, surface acoustic waves, 
and inelastic light scattering will provide other invaluable information inacces- 
sible by the dc magnetotransport. Finally, the real-space imaging ^ of the 
stripes and bubbles would be the most definitive proof of their existence. 

It is a safe bet that we will see a great number of surprises and new devel- 
opments in this area. 
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